Super-horizon perturbations and preheating 
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It has recently been claimed by Bassett et al. that preheating after inflation may affect the amplitude 
of curvature perturbations on large scales, undermining the usual inflationary prediction. We analyze 
the simplest model, and confirm the results of Jedamzik and Sigl and of Ivanov that in linear 
perturbation theory the effect is negligible. However the dominant effect is second-order in the field 
perturbation and we show that this too is negligible, and hence conclude that preheating has no 
significant influence on large-scale perturbations in this model. We briefly discuss the likelihood of 
an effect in other models. 



PACS numbers: 98.80.Cq 



I. INTRODUCTION 
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The standard inflationary paradigm is an extremely 
successful model in explaining observed structures in 
the Universe (see Refs. p],^| for reviews). The inhomo- 
geneities originate from the quantum fluctuations of the 
inflaton held, which on being stretched to large scales be- 
come classical perturbations. The field inhomogencities 
generate a perturbation in the curvature of comoving hy- 
persurfaces, and later on these inhomogeneities are in- 
herited by matter and radiation when the inflaton field 
decays. In the simplest scenario, the curvature pertur- 
bation on scales much larger than the Hubble length is 
constant, and in particular is unchanged during the infla- 
ton decay. This enables a prediction of the present-day 
perturbations which does not depend on the specific cos- 
mological evolution between the late stages of inflation 
and the recent past (say, before nucleosynthesis). 

It has recently been claimed p] that this simple pic- 
ture may be violated if inflation ends with a period of 
preheating, a violent decay of the inflaton particles into 
another field (or even into quanta of the inflaton field it- 
self). Such a phenomenon would completely undermine 
the usual inflationary picture, and indeed the original 
claim was that large-scale perturbations would be ampli- 
fied into the non-linear regime, placing them in conflict 
with observations such as measurements of microwave 
background anisotropies. Given the observational suc- 
cesses of the standard picture, these claims demand at- 
tention. 

In a companion paper [|J , we discuss the general crite- 
ria under which large-scale curvature perturbations can 
vary. As has been known for some time, this is possi- 
ble provided there exist large-scale non-adiabatic pres- 
sure perturbations, as can happen for example in multi- 
field inflation models ■ Under those circumstances 
a significant effect is possible during preheating, though 
there is nothing special about the preheating era in this 



respect and this effect always needs to be considered in 
any multi-component inflation model. 

In this paper we perform an analysis of the simplest 
preheating model, as discussed in Ref. ||. We identify 
two possible sources of variation of the curvature pertur- 
bation. One comes from large-scale isocurvature pertur- 
bations in the preheating field into which the inflaton de- 
cays; we concur with the recent analyses of Jedamzik and 
Sigl H and Ivanov (|] that this effect is negligible due to 
the rapid decay of the background value of the preheat- 
ing field during inflation. However, we also show that 
in fact a different mechanism gives the dominant contri- 
bution, which is second-order in the field perturbations 
coming from short-wavelength fluctuations in the fields. 
Nevertheless, we show too that this effect is completely 
negligible, and hence that preheating in this model has no 
significant effect on large-scale curvature perturbations. 



II. PERTURBATION EVOLUTION 

An adiabatic perturbation is one for which all pertur- 
bations Sx share a common value for Sx/x, where x is 
the time dependence of the background value of x. If 
the Universe is dominated by a single fluid with a defi- 
nite equation of state, or by a single scalar field whose 
perturbations start in the vacuum state, then only adi- 
abatic perturbations can be supported. If there is more 
than one fluid, then the adiabatic condition is a special 
case, but for instance is preserved if a single inflaton field 
subsequently decays into several components. However, 
perturbations in a second field, for instance the one into 
which the inflaton decays during preheating, typically vi- 
olate the adiabatic condition. 

We describe the perturbations via the curvature per- 
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turbation on uniform-density hypersurfaces, denoted CJ3 
In linear theory the evolution of £ is well known, and 
arises from the non-adiabatic part of the pressure per- 
turbations. In any gauge, the pressure perturbation can 
be split into adiabatic and entropic (non-adiabatic) parts, 
by writing 

Sp = c 2 Sp + 6p n M , (1) 
where cj = pj p and the non-adiabatic part is 



<5pnad = pT = p 



5p 5p 



P 



P 



(2) 



The entropy perturbation T, defined in this way, is gauge- 
invariant, and represents the displacement between hy- 
persurfaces of uniform pressure and uniform density. 

On large scales anisotropic stress can be ignored when 
the matter content is entirely in the form of scalar fields, 
and in its absence the non-adiabatic pressure perturba- 
tion determines the variation of £, according to the equa- 
tion §|§ 



dN 



—J- = -iHciT , 



(3) 



where N = In a measures the integrated expansion and 
H is the Hubble parameter. The uniform-density hyper- 
surfaces become ill-defined if the density is not a strictly 
decreasing function along worldlines between hypersur- 
faces of uniform density, and one might worry that this 
undermines the above analysis. However we can equally 
well derive this evolution equation in terms of the den- 
sity perturbation on spatially-flat hypersurfaces, 5p^ = 
— (dp/dN)C,, which remains well-defined. Spatially-flat 
hypersurfaces are automatically separated by a uniform 
integrated expansion on large scales, so the perturbed 
continuity equation in this gauge takes the particularly 
simple form 



dN 



-?>(8pip + 5p.,p) . 



(4) 



From this one finds that Sp^ oc dp/dN for adiabatic per- 
turbations and hence again we recover constant value for 
£. However it is clearly possible for entropy perturbations 
to cause a change in £ on arbitrarily large scales when the 
non-adiabatic pressure perturbation is non-negligible. 



"This is the notation of Bardeen, Steinhardt and Turner |L0j . 
General issues of perturbation description and evolution are 
discussed in a companion paper [Q. The curvature pertur- 
bation of comoving spatial hypersurfaces, usually denoted by 
1Z is practically the same as £ well outside the horizon, 

since the two coincide in the large-scale limit. 



III. PREHEATING 

During inflation, the reheat field into which the infla- 
ton field decays possesses quantum fluctuations on small 
scales just like the inflaton field itself. As these pertur- 
bations are uncorrelated with those in the inflaton field, 
the adiabatic condition will not be satisfied, and hence 
there is a possibility that £ might vary on large scales. 
Only direct calculation can demonstrate whether the ef- 
fect might be significant, and we now compute this effect 
in the simplest preheating model, as analyzed in Ref. [||. 
This is a chaotic inflation model with scalar field poten- 
tial 
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where </> is the inflaton and x the reheat field. Slow- 
roll inflation proceeds with </> > mp\ and gx <C m. The 
effective mass of the \ field is gift and thus will be much 
larger than the Hubble rate, H ~ \J An / '3 rruj) j 'mpi, for 
g ^> m/mp\ ~ 10~ 6 . Throughout, we use the symbol '~' 
to indicate equality within the slow-roll approximation. 

This model gives efficient preheating, since the effective 
mass of x oscillates about zero with large amplitude. In 
most other models of inflation, preheating is less efficient 
or absent, because the mass oscillates about a nonzero 
value and/or has a small amplitude. 

Any variation of £ during preheating will be driven 
by the (non-adiabatic part of) the x field perturbation. 
Our calculation takes place in three steps. The first is 
to compute the perturbations in the x field at the end 
of inflation. The second is to compute how these per- 
turbations are amplified during the preheating epoch by 
the strong resonance. Finally, the main part of the cal- 
culation is to compute the change in £ driven by these x 
perturbations. 



A. The initial quantum fluctuation of the x-Aeld 



Perturbations in the x field obey the wave equation 



S X + 3HS X +['^ + g 2 4> 2 ) S X = 



(6) 



The slow-roll conditions ensure that the x field remains 
in the adiabatic vacuum state for a massive field 



(7) 



where uj z = k 2 /a + g 4> ■ This is a solution provided 



»1, 



(8) 



where m x = g<j) is the effective mass of the x field. 
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The power spectrum of a quantity x, decomposed into 
Fourier components Xk, is denned as 



C. Change in the curvature perturbation on large 
scales 



(9) 



where k = |k| and the average is over ensembles. Hence 
the power spectrum for long-wavelength fluctuations 
(k <C Tn-x) in the \ held simply reduces to the result 
for a massive field in flat space 
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4ir 2 i 



(10) 



where m x is the mass of the field at the required time. 
Physically, this says that at all times the expansion of the 
Universe has a negligible effect on the modes as compared 
to the mass. In particular, at the end of inflation we can 
write 



1 /H, 



X I end 



end 



2tt 



end 



(ii) 



The power spectrum has a spectral index n s = 3. This 
is the extreme limit of the mechanism used to give a blue 
tilt in isocurvature inflation scenarios fllll . 



B. Parametric resonance 

After inflation, the inflaton held <j> oscillates. Strong 
parametric resonance may now occur, amplifying the ini- 
tial quantum fluctuation in \ to become a perturbation 
of the classical field \. The condition for this is 



9 ; 



g 2$2 

~> 1 , 

4m. 2 



(12) 



where $ is the initial amplitude of the 0- field oscillations. 

We model the effect of preheating on the amplitude of 
the x field following Ref. [O] as 



and the Floquet index fi^ is taken as 

^i-ln(l + 2e- 2 ), 

with 



K 2 = 



k \ 2 1 ( k 



18^/g \ fc end 



(13) 



(14) 



(15) 



For strong coupling (q ^> 1), we have k 2 <C 1 for all 
modes outside the Hubble scale after inflation ends (k < 
fcend)- Therefore /ifc ~ ln3/27r 0.17 is only very weakly 
dependent on the wavenumber k. Combining Eqs. (|IT 
and (O) gives 



H, 



end 



2tt 



^enc 



(16) 



In order to quantify the effect parametric growth of the 
X field fluctuations during preheating might have upon 
the standard predictions for the spectrum of density per- 
turbations after inflation, we need to estimate the change 
in the curvature perturbation £ on super-horizon scales 
due to entropy perturbations on large-scales. 

The density and pressure perturbations due to first- 
order perturbations in the inflaton field on large scales 
(i.e. neglecting spatial gradient terms) are of order 
g 2 4> 2 xSx- Not only are the field perturbations Sx 
strongly suppressed on large scales at the end of infla- 
tion [as shown in our Eq. (pd|)] but so is the background 
field x- We can place an upper bound on the size of the 
background field by noting that in order to have slow-roll 
chaotic inflation (dominated by the m 2 (j> 2 /2 potential) 
when any given mode k which we are interested in crossed 
outside the horizon, we require x ^ Tn/g. The large ef- 
fective mass causes this background field to decay, just 
like the super-horizon perturbations, and at the end of in- 
flation we require x ^ m/g(k/k em i) 3 ^ 2 when considering 
preheating in single-field chaotic inflation. Combining 
this with Eq. (11) we find that the spectrum of density 
or pressure perturbations due linear perturbations in the 
X field has an enormous suppression for k <C fc en d : 



^Vxlcnd ^ 



/47T / m \ ( m Pl H CDd 
3 \gmpij \ 2n 



^end 



(17) 



Effectively the density and pressure perturbations have 
no term linear in 5x, because that term is multiplied by 
the background field value which is vanishingly small. 

By contrast the second-order pressure perturbation is 
of order g 2 <fi 2 5x 2 where the power spectrum of Sx 2 is 
given by pi 



fc 3 /- fc -n x (|k'|)n x (|k-k'n 3l , 



2tt 



|k|' 3 |k-k'| 3 



(18) 



We impose the upper limit fc cut ~ k niax to eliminate 
the ultraviolet divergence associated with vacuum state. 
Substituting in for Vsx from Eq. (jl"l|), we can write 



V, 



Sx 2 



8ir / m V 
9 \gmp\J 



H. 



end 



2tt 



^end 



&cnd 



(19) 



Noting that H Qnd ~ m and fc cut ~ fc max — g 1/4 fc C nd, it 
is evident that the second-order effect will dominate over 
the linear term for k < g 1 ^ 2 q 1 ^k end . 

The leading-order contributions to the pressure and 
density perturbations on large scales are thus 
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5p = m 2 (j) 5(j) + 
Sp = —m 2 (f) S(f> 



\g 2 <P 2 5x 2 + \s x 2 , (20) 
-\g 2 tf&x 2 + \5x- (21) 



We stress that we will still only consider first-order per- 
turbations in the metric and total density and pressure, 
but these include terms to second-order in 8\- From 
Eqs. (§), @ and @ we obtain 



<5Pnad = 



-m 2 (j)Sx + (t>g 2 4 ,2 $x 2 



(22) 



where the long-wavelength solutions for vacuum fluctua- 
tions in the <fi field obey the adiabatic condition S4>/4> = 
8(j)/(j). Inserted into Eq. (||), this gives the rate of change 
of C 

Note that the non-adiabatic pressure will diverge pe- 
riodically when <j) = as the comoving or uniform 
density hypersurfaces become ill-defined. Such a phe- 
nomenon was noted in the single-field context by Finelli 
and Brandenberger JT3], who evaded it by instead us- 
ing Mukhanov's variable u = aSip^p which renders well- 
behaved equations. Linear perturbation theory remains 
valid as there are choices of hypersurface, such as the 
spatially-flat hypersurfaces, on which the total pressure 
perturbation remains finite and small. In particular, we 
can calculate the change in the density perturbation due 
to the non-adiabatic part of the pressure perturbation on 
spatially-flat hypersurfaces from Eq. (0), which yields 



A/Onad = -3 / 5p nad Hdt 



(23) 



Even though <5p n ad contains poles whenever = 0, the 
integrated effect remains finite whenever the upper and 
lower limits of the integral are at <$> ^ 0. From this den- 
sity perturbation calculated in the spatially-flat gauge 
one can reconstruct the change in the curvature pertur- 
bation on uniform density hypersurfaces 



AC 



H 



Ap na d 



(24) 



Substituting in our expression for <5p na d we obtain 



AC = ^ 



2m 2 <A 



3H<P 



g 2 cf 2 \5x 2 \Hdt, 



(25) 



where we have averaged over short timescale oscillations 

'2 

of the X - ne ld fluctuations to write 8\ = g 2 4> 2 \&X \- 
To evaluate this we take the usual adiabatic evolution 
for the background <j> field after the end of inflation 



= $ 



sin(mAt) 



mAt ' 

and time-averaged Hubble expansion 



(26) 



H = 



2m 



3 (mAt + 6) 



(27) 



where is an integration constant of order unity. The 
amplitude of the x-held fluctuations also decays propor- 
tional to 1/ At over a half-oscillation from mAt = mr to 
mAt = (n + 1)tt, with the stochastic growth in particle 
number occurring only when <j) = 0. Thus evaluating AC 
over a half-oscillation At = ir/m we can write 



AC = 



2g 2 \S X 2 



3m 2 



1 







(28) 



where x = mAt, s(x) = smx/x, x n = nir and a dash in- 
dicates differentiation with respect to x. The integral is 
dominated by the second term in the bracket which has 
a pole of order 3 when s' — 0. Although s/s' diverges, it 
yields a finite contribution to the integral which can be 
evaluated numerically. For x n 3> 1 the integral is very 
well approximated by , independent of the integra- 

tion constant 0. 

This expression gives us the rate of change of the cur- 
vature perturbation £ due to the pressure of the field 
fluctuations 5\ 2 over each half-oscillation of the inflaton 
field <p. Approximating the sum over several oscillations 
as a smooth integral and using Eq. ( p"3| ) for the growth of 
the x-field fluctuations during preheating (neglecting the 
weak /c-dependence of the Floquet index, fi^, on super- 
horizon scales) we obtain 



c 



ii ad 



16 5 2 |<5x 2 



I end ^2fimAt 



(29) 



The statistics of these second-order fluctuations are 
non-Gaussian, being a ^-distribution. Both the mean 
and the variance of £ n ad are non-vanishing. The mean 
value will not contribute to density fluctuations, but 
rather indicates that the background we are expanding 
around is unstable as energy is systematically drained 
from the inflaton field. We are interested in the vari- 
ance of the curvature perturbation, and in particular the 
change of the curvature perturbation power spectrum on 
super-horizon scales which is negligible if the power spec- 
trum of Cnad on those scales is much less than that of £ 
generated during inflation, the latter being required to 
be of order 10~ 10 to explain the COBE observations. 

To evaluate the power spectrum for £nad we must eval- 
uate the power spectrum of 8\ 2 which is given by sub- 
stituting Vs x , from Eq. (|l6|), into Eq. (18). This gives 



(30) 



where 



/(K,mAt) = ^ f CUt dn' f d6»e 2( ^'+^' )mA V 2 sin6l, 
2 Jo Jo 

(31) 
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k = k/k max as defined in Eq. (|l5|), and is the angle 
between k and k'. Note that at the end of inflation we 
have I(k, 0) = K 3 ut ~ 1, and Ps x 2 ^ This yields 

^Cad - ^5^2 ( mpl ) ( m ) 9 q (fccnd) 

(32) 

One might have thought that the dominant contribu- 
tion to Cnad on large scales would come from 8\ fluctua- 
tions on those scales, and that is indeed the presumption 
of the calculation of Bassett et al. ||. However, in fact 
the integral is initially dominated by kl ~ fc cut , namely 
the shortest scales. The reason for this is the steep slope 
of V$ x ] were it much shallower (spectral index less than 
3/2), then the dominant contribution would come from 
large scales. 

To study the scale dependence of I(k, mAt) and hence 
■PCnad at later times, we can expand /i K - K ' for kk' < 1 as 

2k' cos 9 „, . . 

= M«' + 2 e „ K , 2 k + 0(k 2 ) . (33) 

We can then write the integral in Eq. ([H]) as 

I(k, mAt) = I (mAt) + 0{k, 2 ) , (34) 
where first-order terms, O(k), vanish by symmetry and 

I Q (mAt) = I / KCUt e 4 ^' mA V 2 d K ' . (35) 
2 Jo 

Thus the scale dependence of "P£ nad remains k 3 on large- 
scales for which k«1. 

At late times these integrals become dominated by the 
modes with k' 2 <C (mAt) -1 which are preferentially am- 
plified during preheating. These are longer wavelength 
than fc cu t, but still very short compared to the scales 
which give rise to large scale structure in the present 
Universe. From Eq. ( |l4| ) we have fi K i » /iq — k' 2 /3, for 
k' 2 <C 1, where [Iq = (ln3)/27r, which gives the asymp- 
totic behaviour at late times 

Jo ~ 0.86(mAt)- 3 / 2 e 4w,mAt . (36) 

Thus although the rate of growth of 7\ nad becomes deter- 
mined by the exponential growth of the long-wavelength 
modes, the scale dependence on super-horizon scales re- 
mains proportional to fc 3 for n < (mAt) -1 / 2 . This en- 
sures that there can be no significant change in the cur- 
vature perturbation, £, on very large scales before back- 
reaction on smaller scales becomes important and this 
phase of preheating ends when mAt ~ 100 |l2| ]. 

Numerical evaluation of Eq. ( |32"| ) confirms our analyti- 
cal results, as shown in Fig. 1. For k <§C k max , the spectral 
index remains k 3 during preheating. Observable scales 
have log 10 /c/fc max ~ -20. 
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FIG. 1. The power spectrum of the non-adiabatic curva- 
ture perturbation V( ., shown at four different times: from 
bottom to top mAt = 0, 50, 100 and 150. The parameters 
used were g = 10 -3 , m = 10 _6 mpi and fc cut = fcmax- 

Our result shows that because of the fc 3 spectrum of 5\, 
which leads to a similarly steep spectrum for Cnad, there 
is a negligible effect on the large-scale perturbations be- 
fore the resonance ceases. The suppression of the large- 
scale perturbations in S\, discussed in Refs. |8]|J, means 
that large-scale perturbations in Sx are completely unim- 
portant. However, it turns out that they don't give the 
largest effect, which comes from the short-scale modes 
which dominate the integral for Cnad- Nevertheless, even 
they give a negligible effect, again with a k 3 spectrum. 
Indeed, that result with hindsight can be seen as in- 
evitable; it has long been known |H| that local processes 
conserving energy and momentum cannot generate a tail 
shallower than k 3 (with our spectral index convention) 
to large scales, which is the Fourier equivalent of real- 
izing that in real space there is an upper limit to how 
far energy can be transported. Any mechanism that 
relies on short-scale phenomena, rather than acting on 
pre-existing large-scale perturbations, is doomed to be 
negligible on large scales. 

IV. CONCLUSIONS AND DISCUSSION 

As discussed in detail in a companion paper B , large- 
scale curvature perturbations can vary provided there is 
a significant non-adiabatic pressure perturbation. This is 
always possible in principle if there is more than one field 
or fluid, and since preheating usually involves at least one 
additional field into which the inflaton resonantly decays, 
such variation is in principle possible. 

In this paper we have focussed on the simplest preheat- 
ing model, as discussed in Ref. 0. We have identified 
the non-adiabatic pressure, and shown that the domi- 
nant effect comes from second-order perturbations in the 
preheating field. Further, the effect is dominated by per- 
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turbations on short scales, rather than from the resonant 
amplification of non-adiabatic perturbations on the large 
astrophysical scales. Nevertheless, we have shown that 
the contribution has a k 3 spectrum to large scales, ren- 
dering it totally negligible on scales relevant for structure 
formation in our present Universe by the time backreac- 
tion ends the resonance. Amongst models of inflation 
involving a single-component inflaton field, this model 
gives the most preheating, and so this negative conclu- 
sion will apply to all such models. 

Recently Bassett et al. Q have suggested large effects 
might be possible in more complicated models. They 
consider two types of model. In one kind, inflation takes 
place along a steep-sided valley, which lies mainly along 
the direction of a field cf> but with a small component 
along another direction x- I n this case, one can simply 
define the inflaton to be the field evolving along the val- 
ley floor, and the second heavy field lies orthogonal to it. 
Taking that view, there is no reason to expect the pre- 
heating of the heavy field to give rise to a bigger effect 
than in the simpler model considered in this paper. 

In the second kind of model, the reheat field is 
light during inflation, and this corresponds to a two- 
component inflaton field. As has long been known, there 
can indeed be a large variation of £ in this case, which can 
continue until a thermalized radiation-dominated uni- 
verse has been established. Indeed, in models where one 
of the fields survives to the present Universe (for exam- 
ple becoming the cold dark matter), variation in £ can 
continue right to the present. This variation is due to 
the presence on large scales of classical perturbations in 
both fields (properly thought of as a multi-component 
inflaton field) generated during inflation, and the effect 
of these must always be considered in a multi-component 
inflation model, with or without preheating. 
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